Is the Broad-Line Region Clumped or Smooth? 
Constraints from the Ha Profile in NGC 4395, the Least 
Luminous Seyfert 1 Galaxy 
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ABSTRACT 

The origin and configuration of the gas which emits broad lines in Type I active galactic 
nuclei is not established yet. The lack of small-scale structure in the broad emission- line profiles 
is consistent with a smooth gas flow, or a clumped flow with many small clouds. An attractive 
possibility for the origin of many small clouds is the atmospheres of bloated stars, an origin which 
also provides a natural mechanism for the cloud confinement. Earlier studies of the broad- line 
profiles have already put strong lower limits on the minimum number of such stars, but these 
limits are sensitive to the assumed width of the lines produced by each cloud. Here we revisit this 
problem using high-resolution Keck spectra of the Ha line in NGC 4395, which has the smallest 
known broad- line region (~ 10 cm). Only a handful of the required bloated stars (each having 
7'* ~ 10 14 cm) could fit into the broad-line region of NGC 4395, yet the observed smoothness of 
the Ha line implies a lower limit of ~ 10 4 — 10 5 on the number of discrete clouds. This rules out 
conclusively the bloated-stars scenario, regardless of any plausible line-broadening mechanisms. 
The upper limit on the size of the clouds is ~ 10 12 cm, which is comparable to the size implied 
by photoionization models. This strongly suggests that gas in the broad- line region is structured 
as a smooth rather than a clumped flow, most likely in a rotationally dominated thick disk-like 
configuration. However, it remains to be clarified why such a smooth, gravity-dominated flow 
generates double-peaked emission lines only in a small fraction of active galactic nuclei. 



Subject headings: galaxies: active 
NGC 4395 

Introduction 
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Type I active galactic nuclei (AGNs) are defined 
by the presence of permitted lines with a kine- 
matically distinct broad component, usually with 
full width at half-maximum intensity (FWHM) 
> 1000 km s _1 . The broad lines are emitted 
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by gas in the broad-line region (BLR) which is 
being photoionized by the central compact con- 
tinuum source. Detailed photoionization calcula- 
tions yield rather tight constraints on the density, 
temperature, ionization level, and chemical abun- 
dances of the photoionized gas (e.g., Ferland et al. 
1998; Hamann & Ferland 1999), while reverber- 
ation mapping techniques constrain the size and 
velocity field of the gas (e.g., Peterson 2001). How- 
ever, the structure and origin of the photoionized 
gas in the BLR remains a largely unsolved prob- 
lem. 

1.1. Proposed Models for the BLR 

The gas in the BLR may be structured in small 
clumps (hereafter "clouds"). A particularly ap- 
pealing model is that of bloated stars, where the 
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clouds are identified with the surface layers of su- 
pergiant stars, which are further bloated by the 
ionizing radiation field (Edwards 1980; Scoville & 
Norman 1988; Pcnston 1988; Kazanas 1989; Tout 
ct al. 1989; Alexander & Netzer 1994, 1997). This 
model is attractive because it relies on a relatively 
small modification of known objects (supcrgiants), 
and it also provides a natural mechanism for the 
cloud confinement, although a self-consistent BLR 
structure within this model appears problematic 
(Begelman & Sikora 1992; Krolik 1999). 

Alternative suggestions for the origin of BLR 
clouds include gas streams produced by tidally 
disrupted stars (Hills 1975; Roos 1992), or by 
star-disk collisions (Zurek et al. 1994; Vilko- 
viskij & Czerny 2002), or clumps in a gravita- 
tionally unstable outer disk (Collin & Hure 2001). 
Other scenarios invoke density inhomogeneities 
produced by wind interaction with nearby obsta- 
cles (Perry & Dyson 1985), such as stellar atmo- 
spheres (Torricelli-Ciamponi & Pietrini 2002), su- 
pernova remnants (Pittard et al. 2003), or the 
the accretion-disk surface (Cassidy & Raine 1996). 
Alternatively, it was suggested that the clouds are 
density enhancements produced by shocks in an 
accreting gas (Fromerth & Melia 2001). The phys- 
ical viability of some of the proposed mechanisms 
remains unclear (e.g., Mathews 1986; Mathews & 
Doane 1990) 

The other option for the BLR structure is a 
smooth continuous flow, rather than a clumped 
flow. For example, this flow may be identified with 
the illuminated and photoionized accretion-disk 
surface (e.g., Dumont & Collin-Souffrin 1990a, b). 
Alternatively, it could originate in a wind driven 
off an accretion disk by radiation or magnetic pres- 
sure gradients (Shlosman et al. 1985; Emmcring 
ct al. 1992; Konigl & Kartje 1994; Murray et 
al. 1995; Chiang & Murray 1996; Bottorff ct al. 
1997). 

1.2. Earlier Constraints 

The minimum size of the photoionized clouds 
can be estimated from simple Stromgren depth ar- 
guments. Equating the H ionization rate to the re- 
combination rate per unit area yields the column 
density of the photoionized H II layer, where most 
of the line emission is produced. This argument 
yields S; on 10 23 C/ cm~ 2 , where U = nj/n e is 
the ionization parameter, and n 7 and n e are (re- 



spectively) the ionizing photon and the electron 
number densities. Typical BLR parameters are 
U ~ 0.1 and n e 10 10 cm -3 , implying a typical 
thickness of the H II layer of d « 10 12 cm (al- 
though there may be a large range of d values in 
the BLR; e.g., Baldwin et al. 1995). The equiv- 
alent widths of the H recombination lines imply 
that the BLR gas covers f2BLR./47r rs 0.1 — 0.2 
of the ionizing continuum source (e.g., Ferland et 
al. 1998). The photoionized gas could be made 
of clumps which are as small as d, or it may re- 
side in a thin photoionized surface layer over more 
extended bodies. The latter option is supported 
by the strength of some of the low-ionization lines 
(e.g., lines from Fe II, Ca II, Na I) which originate 
in deeper layers at £ > Ei on , where H is only par- 
tially ionized (e.g., Ferland & Persson 1989; Bald- 
win et al. 2003). How can one then determine if 
the BLR is clumped or smooth? 

The number of clouds in the BLR, n Cl and their 
size, r Cl can be constrained based on the smooth- 
ness of the emission-line profiles. The level of 
fluctuations in the emission-line profiles is pro- 
portional to 1/ y/nZ, assuming the cloud velocities 
are randomly distributed within the line profile 
(in proportion to the line flux). The solid area 
subtended by all the BLR clouds is Oblr/47t = 
n c 7rr 2 /47ri?g LR , where -Rblr is the radius of the 
BLR; since Qblr/47t w 0.1 - 0.2, we find that 
r c ~ -Rblr/v^c- High-resolution spectroscopy at 
a high signal-to- noise ratio (S/N) can place tight 
constraints on the level of fluctuations, which then 
translates into a large lower limit on n c , and thus 
a tight upper limit on r Cl which may exclude some 
clumped BLR models. 

The first application of this method was made 
by Capriotti et al. (1981), who deduced n c £ 10 4 
in a sample of Seyfert galaxies. Atwood et al. 
(1982) improved the limit to n c £ 5 x 10 4 for 
Mrk 509, and Laor et al. (1994) obtained n c > 10 4 
for a sample of bright quasars, using Hubble Space 
Telescope spectroscopy of the C IV and Ly a lines. 
Significantly improved limits of n c ^ 3 x 10 6 and 
3 x 10 7 were obtained for Mrk 335 (Arav et al. 
1997) and NGC 4151 (Arav ct al. 1998) using 
high-quality Keck spectra of the Ha emission line. 
Finally, Dietrich et al. (1999) obtained the highest 
limit of n c £ 10 s in the luminous quasar 3C 273. 

Since i? B LR oc L~ p ° t 6 (Kaspi ct al. 2005), while 
r c f« -Rblr/V^ci tighter constraints on r c can be 
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obtained in lower-luminosity objects. In 3C 273, 
jaL„(5100A) w 6 x 10 45 erg s"\ and the low- 
est luminosity object where n c was constrained is 
NGC 4151, where j/L y (5100A) = 7 x 10 42 erg s" 1 
(Kaspi et al. 2000). 

Here we report on high-quality Keck spec- 
troscopy of the Ha line in NGC 4395, which is 
a factor of 1000 fainter than NGC 4151, having 
^(5100A) = 5.9 x 10 39 erg s" 1 (Peterson et al. 
2005). This is the least luminous known Seyfert 1 
nucleus (Filippenko & Sargent 1989; Filippenko, 
Ho, & Sargent 1993), and despite its extreme low 
luminosity, it appears to have a rather normal 
broad emission-line spectrum (apparently with a 
lower than solar metal abundance, Kraemer et al. 
1999). The low luminosity of NGC 4395 implies 
an expected -Rblr of just ~ (0.5 — 1.7) x 10 14 cm, 
based on the most recent expression for the mean 
Balmer-line time lag versus ^L„(5100A) in Kaspi 
et al. (2005; see their Table 3). Bloated stars are 
predicted to have a size of ~ 10 14 cm (Alexander 
& Netzer 1994), and thus at most a handful of 
such stars can be accommodated in the BLR of 
NGC 4395. This will inevitably produce large- 
amplitude features in the line profile. 

The earlier studies of Arav et al. (1997, 1998) 
and Dietrich et al. (1999) ruled out the bloated- 
stars idea, assuming the width of the lines pro- 
duced by the individual clouds are not much above 
the thermal width. However, significant nonther- 
mal broadening through radiation pressure abla- 
tion of the BLR clouds, or through magnetic pres- 
sure effects, is likely and possibly inevitable. This 
broadening reduces the lower limit on n c (e.g., 
Bottorff & Ferland 2000; Bottorff et al. 2000), 
and it may leave the bloated-stars idea as a viable 
option. As we show below, the bloated-stars idea 
can be conclusively ruled out in the compact BLR 
of NGC 4395, regardless of the velocity width of 
the lines produced by the individual clouds. The 
upper limit obtained is comparable to d, strongly 
suggesting that the broad lines are produced by a 
smooth gas flow. In §2 of this paper we describe 
the observations and method of analysis, and the 
results are discussed in §3. 

2. Observations and Analysis 

The spectra we use here were obtained in three 
observations of NGC 4395 with the Keck-I and 



Keck-II 10 m telescopes. Two of the observations 
were made with the High Resolution Echelle Spec- 
trometer (HIRES; Vogt et al. 1994); the spectral 
resolution was A/AA = 38,000, or 0.173 A near 
Ha, sampled at 0.0473 A per pixel. The third ob- 
servation was made with the Echellette Spectro- 
graph and Imager (ESI; Shcinis et al. 2002); the 
spectral resolution was A/AA = 8,000, or 0.820 A 
near Ha, sampled at 0.259 A per pixel. The first 
HIRES observation (hereafter HIRES1), described 
by Filippenko & Ho (2003), was a 20 + 40 minute 
exposure obtained on 1994 April 14 UT, which 
yielded a S /N ranging from 20 per pixel at the con- 
tinuum near Ha to 300 at the line peak. The sec- 
ond HIRES observation (hereafter HIRES2) was a 
30 + 60 minute exposure taken on 2002 February 
4 UT, yielding a S/N ranging from 30 to 300, as 
above. The ESI observation was a 20 + 20 minute 
exposure made on 2002 December 2 UT, and the 
S/N ranges from 50 to 400, as above. 

All of the spectra were reduced using the MA- 
KEE reduction software 1 . Cosmic-ray hits were 
manually removed from the spectra by comparing 
the two exposures available for each of the spec- 
tra. We did not use standard cosmic-ray removal 
tasks as they sometimes leave small residual fea- 
tures, which will affect the fluctuation analysis. 
We used a simple summation for the extractions 
rather than an optimal extraction in order to pre- 
serve the correct emission-line profiles, since an 
optimal extraction might produce incorrect results 
for the narrow-line profiles, and it could poten- 
tially produce small distortions in the broad lines 
as well. 

The rest-frame wavelength scale was converted 
to a velocity scale with respect to the rest wave- 
length of Ha, 6562.8 A. We applied a slight off- 
set to the velocity scale in each spectrum, such 
that the narrow Ha peak falls at v = km s~ 4 in 
all cases. We next subtracted a constant contin- 
uum level from each spectrum, in order to get the 
net line flux. Absolute flux calibration is available 
only for the ESI spectrum, and an arbitrary flux 
scale is used for the two HIRES spectra. The lack 
of an absolute flux calibration does not affect our 
analysis as we are only interested in the Ha profile 
shape here. 

The minimum number of clouds is estimated 
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as follows. We first find for each observed profile, 
fx hs , the smoothest acceptable fit, — that is, 
a fit which yields Xr ~ 1 f° r the given flux mea- 
surement errors, /^ rr . We next produce a Monte 
Carlo (MC) realization of the observed profile by 
adding the emission of n c lines, each one having a 
Gaussian line shape with a velocity dispersion a. 
Each line is further convolved with the instrumen- 
tal resolution, approximated as a Gaussian with 
FWHM = 8 km s" 1 for HIRES, and 38 km s" 1 
for ESI. We then measure % 2 01 the MC simulated 
profile, fx™, with respect to the smooth-fit profile, 
as a function of n c and a. For a given a, we find 
the minimum n c required to obtain x 2 ~ 1- This 
procedure provides us with an estimate of the min- 
imum n c , as a function of a, required so that the 
level of fluctuations produced by the finite num- 
ber of clouds does not exceed the observed level of 
fluctuations (produced by the flux measurement 
errors). 

To obtain the smooth model fit, each line is 
fit with a spline function, where we iterate man- 
ually over the number and positions of the points 
selected for the spline fit, until the fit yields 
a reduced Xr ~ 1 over the velocity range of 
—2500 km s _1 to 2500 km s _1 (or smaller, as set 
by the available spectra). The error spectrum, 
./^ rr , obtained together with ,/° bs in the MAKEE 
pipeline reduction, is used for calculating the x 2 
of the spline fit. The broad Ha line is blended 
with the narrow Ha component and with the nar- 
row [N n] A6548 and [N II] A6583 doublet lines. To 
make sure that the spline fit is made only in re- 
gions dominated by the broad-component emis- 
sion, we add the additional constraint that the 
residual flux in the two narrow [Nil] line pro- 
files should match. 2 The [Nil] profile was also 
used as a guide for the expected shape of the 
narrow Ha component, and the broad Ha spline 
fit was constructed to yield a narrow Ha pro- 
file similar (but not identical) to the [Nil] pro- 
file. This procedure allowed us to determine the 
velocity ranges affected by narrow-line contamina- 
tion, which we find are —900 km s _1 < v < —550 
km s"\ -210 km s _1 < v < 230 km s" 1 , and 700 
km s _1 < v < 1200 km s^ 1 . We exclude these 
velocity ranges in the % 2 calculation. 

2 The measured flux ratio of [Nil] A6583/[Nn] A6548 in the 
spectra is ~ 3, which agrees with the theoretically predicted 
value 



Figure 1 displays the HIRES 1, HIRES2, and 
ESI Ha profiles, the smooth fits, and the residu- 
als from the fit. Achieving x 2 ~ 1 required about 
14 spline points per spectrum. These were spaced 
about 150 km s _1 apart close to the line core 
(between the [Nil] lines), and about 250 km s _1 
or more apart in the line wings. The x 2 & n d the 
number of degrees of freedom obtained through 
the spline fit are 432/400 for ESI, 1118/1446 for 
HIRES1, and 688/812 for HIRES2. Thus, the fluc- 
tuations on scales smaller than the spline smooth- 
ing length are consistent with the level expected 
from /^ rr . There is no evidence for significant in- 
trinsic profile structure beyond pure noise. 3 There 
are significant, though weak (amplitude <~ 1%), 
features in the broad-line profiles close to the core 
on scales ^ 150 km s _1 , as indicated by small 
and non-monotonic changes in the smooth-fit pro- 
file curvature. At the wings the smooth fit ap- 
pears featureless; however, since the S/N drops, 
such broad weak features would not be detectable 
there. 

For the MC simulation, each one of the n c 
clouds is assumed to emit a line with a Gaussian 
profile with the same total flux and velocity dis- 
persion a. The clouds differ only in their centroid 
velocity v. We next need the velocity probability 
distribution function, P(v), for the MC simula- 
tion. This P(v) should give fl im -> in the 
limit of n c — > oo. As an estimate for P(v) we use 
the inverse function of 

/V 
fvdv 
-oo 

(C is set such that = 1, as required for the in- 
tegral of a distribution function) , where /„ is /f ' 
converted to a velocity scale. The use of an in- 
verse function is a general procedure allowing one 
to transform a uniform random distribution from 
to 1 (generated by any standard random-number 
generator) into any given normalized distribution 
(/„ here; see also Press ct al. 1992, Chapter 7.2). 
The MC profile, fx™ 1 , is generated by adding the 
contribution of the n c individual profiles, each one 
at its respective v. The goodness of fit for each 
MC realization is measured by calculating the x 2 
of fx™ 1 , with respect to /f *, using /| rr for the flux 

3 This also rules out the presence of detectable stellar ab- 
sorption features from the underlying host galaxy light. 
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errors. Since each MC realization yields a some- 
what different value of x 2 , we performed 20 dif- 
ferent MC simulations for each set of n c and a 
values, and obtained the average value \r av an( ^ 
Ax 2 ., the dispersion of xl in the 20 simulations, 
which is used for verifying that Xr,av 1S accurate 
enough (i.e., Ax 2 /Xrav ^ !)■ A model was de- 
fined as acceptable if x 2 av = 1 ± A% 2 . 

The observed line profile is a convolution of the 
cloud velocity distribution function, P(v), and the 
line profile produced by each cloud. Above we 
have made the approximation that P(v) is propor- 
tional to the observed profile, and the accuracy of 
this approximation drops as a increases. Specifi- 
cally, we find that for b > 50 km s _1 (b = \plu) 
the above procedure does not yield /| lm — > / A * 
for n c — > oo. The simulated line peak becomes 
broader than the observed peak due to the effect of 
the convolution. Another difficulty arises from the 
fact that the spline fit we use follows profile fea- 
tures on scales £ 150 km s _1 , and thus features 
on this scales are filtered out of the observed resid- 
uals spectrum (f x hs — fx*)- On the other hand, 
a MC simulation of clouds with a given a pro- 
duces spectral features on a scale of a. Thus, for 
large a, fx" 1 — /f* shows large-scale features, while 
such features do not exist in the observed residuals 
spectrum. These features can produce significant 
X 2 even for large values of n c . To overcome these 
two problems we added the following iteration. We 
took /| lm obtained for a specific simulation, and 
fit it with a spline (using the same velocities used 
in the observed line fit), in order to get a revised 
estimate of /jj*. We now use the revised /f*, to 
measure the Xr of /f lm with respect to the revised 

ffit 

Figure 2 demonstrates the effect of n c and a 
on the level of fluctuations in the MC simula- 
tions, compared to the observed residuals (shown 
for the HIRES2 spectrum). As expected, the sim- 
ulated profile becomes smoother with increasing 
n c and a. In particular, for b = 10 km s _1 , 
which corresponds to pure thermal broadening of 
T w 10 4 K gas, the MC fluctuations are compa- 
rable to the observed ones for n c = 10 5 , while for 
b = 100 km s _1 , n c = 3 x 10 3 is enough to pro- 
duce fluctuations which are about comparable to 
the observed ones. A more quantitative measure 
of the fluctuations level is obtained by calculating 
the x 2 °f the fluctuations, as mentioned above. 



Figure 3 presents our derived limits on the min- 
imum n c , as a function of b, required to yield 
Xl < 1 for the HIRES1/2 spectra and for the ESI 
spectrum. The strongest constraints are provided 
by the ESI spectrum, where the broad Ha line 
happens to have the largest strength relative to 
the narrow Ha line. For b = 10 km s _1 the mini- 
mum number of clouds required is 1.35 x 10 5 ; this 
drops to 2.4 x 10 4 for b = 100 km s"\ and to 
3900 for b = 200 km s" 1 . The solid lines in the 
plot follow the n c cx 6 _1 relation, expected from 
simple fluctuation arguments. 4 The ESI limits at 
b < 40 km s _1 are nearly independent of b be- 
cause the convolution with the instrumental line 
profile (b = 22.7 km s _1 ) leads to a nearly con- 
stant line profile. The upper limits on n c fall below 
the n c cx b^ 1 relation at b > 50 km s _1 in all spec- 
tra due to the profile smoothing procedure, which 
filters out large-scale profile structure. This effect 
is inevitable given the lack of a priori knowledge 
about the intrinsic line profile. We note, however, 
that when b > 200 km s" 1 the simulated profile 
becomes distinctly different from the observed one. 
Our limits are also subject to the assumption that 
all clouds produce similar line profiles. Otherwise, 
any observed profile can be decomposed to any 
number of different profile components, but such 
an arbitrary decomposition appears implausible. 

We note in passing that the strongest constraint 
on n c (that is, the largest contribution to Xr) 
comes from the observed fluctuations in the line 
core, rather than the wings (cf. Dietrich et al. 
1999; Ferland 2004). This arises from the fact 
that the line flux density in the far wings be- 
comes smaller than the continuum flux density 
(here f[ inc > /^ ont for -900 km s _1 < v < 
600 km s" 1 in HIRES1, and for -1100 km s _1 < 
v < 1100 km s- 1 in HIRES2 and ESI). Thus, the 
far-wings noise is dominated by continuum pho- 
tons rather than line photons. Since the contin- 
uum level is constant, the noise level is constant, 
and thus the observed S/N is cx f\ mc (put in other 
terms, the S/N in the far wings is background lim- 
ited, rather than source limited). In contrast, the 

4 The number of clouds contributing to each velocity bin in 
the spectrum is <x n c and also cx b, i.e. it is cx b X n c . The 
number of clouds in each velocity bin sets the fluctuation 
level in the spectrum. Thus, at a fixed level of fluctuations, 
for example the one which corresponds to Xr = b onc 
expects n c cx b^ 1 . 
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MC simulation fluctuations are anywhere propor- 
tional to the square root of the number of clouds 

contributing to a given bin, and so are oc \J /| inc . 

Thus, the ratio of available S/N to the MC simula- 
tion fluctuations, which indicates the detectability 
of the predicted fluctuations, is oc y/fx; i.e. it ap- 
proaches zero at the wings. 

3. Discussion 

The earlier studies of Arav et al. (1997, 1998) 
and Dietrich et al. (1999) have already ruled out 
the bloated-stars scenario based on the lack of the 
expected level of fluctuations in the line profiles. 
However, the fluctuation level is reduced if the line 
width generated by each cloud is much larger than 
the thermal width (c s w 10 km s _1 ), leaving the 
possibility of a clumped BLR if a » c s . Such a 
large a may be produced by, for example, radia- 
tive ablation effects (e.g., Scoville & Norman 1995; 
Pier & Voit 1995). 

The predicted level of fluctuations is cx 
Since n c ss (-RblrAc) 2 , the largest fluctuations 
are expected in objects with the smallest -Rblr — 
that is, in the lowest-luminosity objects. Here we 
report on analysis of Keck observations of the Ha 
profile in NGC 4395, which is a factor of <~ 10 3 less 
luminous than the lowest-luminosity object pre- 
viously analyzed (NGC 4151, Arav et al. 1998), 
and is the lowest-luminosity type I AGN known 
by a significant factor (compare with Barth et al. 
2004; Greene & Ho 2004). Peterson et al. (2005) 
recently measured a size of 1.0 x 10 14 cm for the 
C IV emitting region in NGC 4395. The Balmer 
line emitting region is measured to be about twice 
as large as the C IV emitting region in the few ob- 
jects were both lines were measured (Korista et al. 
1995; Onken & Peterson 2002), and comparison of 
the Ha and C IV line widths suggests it could be 
as much as 5 times larger here, so one may expect 
-Rblr ~ (2 - 5) x 10 14 cm for the Ha line. Thus, 
if the bloated stars have r c w 10 14 cm (Alexander 
& Netzer 1994), then only a handful of them can 
fit in the BLR of NGC 4395, and this will produce 
large profile features, essentially regardless of the 
line width generated by each cloud. Such a model 
is clearly ruled out by the observed smoothness 
of the Ha line, which implies n c £ 10 4 — 10 5 for 
few 10- 100 km s" 1 . 

The presence of r c w 10 14 cm stars at such a 



small distance from the black hole can also be 
ruled out based on tidal disruption arguments, 
as noted by Alexander & Netzer (1997). A star 
with a mass M* and radius r* will be tidally dis- 
rupted once it crosses inward of r t idai ~ r * x 
(Mbh/M,) 1 / 3 . In NGC 4395, M BH = 3.6xlO 5 M 
(Peterson et al. 2005), implying rtidai/'** ~ 70, 
and thus bloated stars, which presumably have 
M* w M Q 5 , will be disrupted already at r t idai ~ 
7 x 10 15 cm. Stars which survive down to -Rblr 
without being tidally disrupted need to be smaller 
than (3 - 7) x 10 12 cm, or (40 - 100) R & ; they 
can at most be red giants, or blue supergiants. 
A minimum population of ~ 5000 such stars is 
required in NGC 4395 to produce the necessary 
BLR covering factor (using Rblr/v* > 70). Such 
a population will have a bolometric luminosity of 
- 10 39 - 10 40 erg s- 1 (e.g., Table 3.14 of Bin- 
ney & Merrifield 1998), and its /-band magnitude 
will be comparable to (or larger than) the unre- 
solved point-source magnitude measured by Filip- 
penko & Ho (2003). In any case, physical colli- 
sions will destroy such an extremely dense stellar 
system (> 10 15 M Q pc~ 3 ) in a few dynamical times 
(a few weeks). Thus, a stellar origin for the BLR 
clouds in NGC 4395 is ruled out both on theo- 
retical grounds, and directly through the profile 
fluctuation analysis. 

The tidal disruption argument can also be used 
to place a lower limit on the density of any self- 
gravitating cloud at the BLR, regardless of its size. 
A cloud of density p c and mass M c = 4np c rl/3 
can avoid tidal disruption at the BLR if p c > 
3/47tMbh/-Rblr- This gives a minimum number 
density, N c > 10 18 enr 3 in NGC 4395, about 10 8 
times larger than the density indicated by pho- 
toionization models (Kraemer et al. 1999). Thus, 
any clumped emission-line gas in the BLR will be 
subject to the strong shearing effect of the black 
hole gravity, which will stretch clouds by Ar c w r c 
within one dynamical time. 

Our lower limit on n c of ~ 10 4 — 10 5 implies 
r c ~ (1 — 3) x 10 12 cm, which is comparable to 
d, the minimum size implied by photoionization 
models. Since wc do not have a robust indica- 
tion here that r c <C d, the BLR gas may still be 

Significantly more massive stars are unlikely due to their 
short lifetime, while compact stars are unlikely due to their 
high surface gravity. 
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composed of clumps with the minimum acceptable 
size. However, it is not clear what mechanism can 
create, confine, and maintain such small individ- 
ual clouds in a steady state (e.g., Done & Krolik 
1996; Krolik 1999). Thus, the most plausible sce- 
nario appears to be a smooth, or quasi-smooth, 
gas flow. Reverberation mapping indicates that 
the BLR Balmer line velocity is dominated by 
gravity, and a pure inflow or outflow is excluded 
(e.g., Peterson 2001). Thus, the only remaining 
velocity field appears to be a rotation-dominated 
flow, most likely in a geometrically thick config- 
uration (to subtend the required covering factor). 
Direct evidence for an outflowing warm ionized gas 
is present in UV (Crenshaw et al. 2004) and X-ray 
observations (Shih et al. 2003), and this outflow 
may originate in the inner denser BLR component. 

A rotation-dominated gas flow should produce 
double-peaked lines. Such lines are indeed ob- 
served in AGNs with very broad lines (e.g., Era- 
cleous & Halpcrn 2003; Strateva et al. 2003), but 
the majority of AGNs show single-peaked lines. 
Chiang & Murray (1996) suggested that a suffi- 
ciently large radial velocity gradient can modify 
the line emissivity and produce a single-peaked 
emission line. However, lower optical depth lines, 
such as semi-forbidden lines and the higher H- 
series lines, should still show double-peaked struc- 
ture, unlike what is generally observed. 

We note in passing that indirect indications for 
strongly blended "hidden" double peaks are pro- 
vided by the rotation of the polarization angle 
across the Ha line profile observed in some objects 
(Smith et al. 2005), and by the distinct variabil- 
ity characteristics of the red and blue line wings 
observed sometimes (e.g., Veilleux & Zheng 1991; 
Wanders & Peterson 1996). Both effects were in- 
terpreted as evidence that the Balmer lines orig- 
inate in a rotating disk. The lack of a double- 
peak structure could be induced by a large ve- 
locity broadening mechanism in the BLR, which 
makes the double peaks blend into a single peak, 
although the double-peak structure would still 
be discernible through the different polarization 
and variability characteristics of the two blended 
peaks. The two peaks appear to separate out in 
objects where the Balmer lines are broad enough 
(e.g., Strateva et al. 2003). 

To summarize, although the BLR gas is likely 
to reside in a smooth flow, the origin of the gas, 



the structure of the flow, and its dynamics remain 
to be clarified. 
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Fig. 1. — The three Keck spectra of the Ha line used for the analysis. The upper part of each subpanel 
shows the observed spectrum (histogram), and the spline fit (smooth line). The lower part of each subpanel 
shows the observed/spline-fit flux ratio; the hatched regions are affected by the narrow emission lines and 
excluded from the fluctuation analysis. Note that the ESI spectrum has the highest S/N ratio, although it 
has a lower spectral resolution compared to the two HIRES spectra. The two HIRES spectra are not flux 
calibrated; they were arbitrarily offset so that all spectra have a matching peak flux density. 
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Fig. 2. — A comparison of the observed versus simulated Ha profile fluctuations. The dotted lines 
show the observed/spline-fit flux ratio for the HIRES2 spectrum. The solid lines shows the Monte Carlo 
simulation/spline- fit flux ratio. Each row shows simulations with the same number of clouds n c , and each 
column shows simulations with the same Gaussian b parameter. Note that at b = 10 km s _1 (a pure thermal 
width), n c ^ 10 5 is required to suppress the fluctuations to the observed level; at b = 30 km s _1 , n c > 10 4 
is required; and at b = 100 km s , 10 3 is sufficient. 
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Fig. 3. — The minimum number of clouds required to suppress the fluctuations to the observed level, as 
a function of the width of the line emitted by each cloud. The tightest constraints are obtained from the 
ESI spectrum, where the broad Ha line has the highest S/N ratio (see Fig. 1). The solid lines follow 
n c oc the expected relation from simple fluctuation arguments. The ESI limits at b < 40 km s _1 are 
nearly independent of b because the convolution with the instrumental line profile leads to a nearly constant 
effective b. At b > 50 km s _1 , the simulation constraints start falling below the —1 slope relation. This is 
a result of the spline smoothing which filters out fluctuations on scales £ 100 km s _1 . 
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